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Abstract 

We study the subcritical bubble formation near the phase space domain wall. We 
take into account that the phase of the scalar field can vary using complex U(l) 
symmetric field and a phenomenological potential with cubic term responsible 
to symmetry breaking. We show that the presence of the domain wall induces 
subcritical bubbles so that their formation rate near the wall is considerably 
larger than far of it. The allowed deviations of the phases of new bubbles are so 
large that they prevent the system from induced nucleation. 
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The study of cosmological phase transitions have received considerable attention 
in literature since the possibility of electroweak baryogenesis was introduced Q. In 
that context the main effort has been performed considering the critical configuration 
(bubble) formation which are stable, expanding broken phase domains by their own. 
Less attention has been paid to the vast subcritical bubble formation during the period 
of the metastable state, because they have been supposed to have no effect to the 
phase transition itself. Also it has become clear 0, |3], Q that for physical Higgs mass, 
heavier that 60 GeV, the electroweak phase transition is only weak. In the context of 
electroweak phase transition, subcritical bubbles were discussed first time by |J and it 
has been also shown || that phase equilibrium can be reached at weak enough phase 
transitions. Kinetics of subcritical bubbles has been studied in where phase mixing 
above the critical temperature was investigated. More recently numerical simulations 
of phase mixing in a 2+1 -dimensional model || and in 3+1 -dimensional model || 
have been performed. 

When the electroweak phase transition is weak enough and the formation rate 
of broken phase configurations is large enough, there are large number of subcritical 
bubbles present at and below the critical temperature. It can happen that a collec- 
tion of subcritical bubbles form a region with a size of a critical bubble before they 
shrink away [|K]]. Therefore the phase transition may be triggered by clustering the 
broken phase regions. This mechanism of producing phase transition is called induced 
nucleation. 

In Ref. |1(J it was demonstrated that in a simple model induced nucleation is 
indeed possible for a large range of parameter values of the scalar potential. The 
growth of a (spherical) region of broken phase was approximated to happen layer by 
layer by nucleating subcritical bubbles. As noted there, however, the probability that 
the region of broken phase grows to the size of a critical bubble strongly depends 
on the relative phases of the subcritical bubbles. This effect was accounted in ||10| 



by introducing a parameter which describes the probability that the phase of the 
nucleated subcritical bubble is correlated with the phase of the pre-existing bubble. 
Its origin, however, was not further analyzed there. 

In the present paper we analyze the subcritical nucleation near phase boundary. 
In particular, we concentrate to the question, how phases of newly formed bubble and 
pre-existing broken phase domain correlate. We assume that the relative velocities 
between the bubbles are negligible. In opposite case, the relative phase of colliding 
bubbles would be washed out by thermal processes with a rate which is presumably 
much larger than in the case of bubbles in rest, leading effectively to the case where 
the phases of the bubbles are correlated. 

One can also consider the effect of the distance u of the subcritical bubble from 
the pre-existing one. If u is small, the newly formed bubble overlaps significantly with 



the pre-existing one, causing the spatial phases to be correlated. As u increases, the 
spatial phase difference becomes more probable. 

We work with a 3 + 1 -dimensional phenomenological model with a complex scalar 
$ having the 3 -dimensional action 

S 3 [$] = J d 3 xl\V<S>\ 2 + V($)], (f) 

where the phenomenological potential is given by 

V[$] = m(T) 2 |$| 2 - ^«T|$| 3 + A|$| 4 . (2) 
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Generally, the non-zero maximum V- and minimum v + of the potential (below the 
critical temperature T c ) are given by 



where 

The parameter A is less than unity at temperatures below T c . 

Generally there exist subcritical bubbles with various sizes and forms. We assume 
that a typical subcritical bubble has a Gaussian, spherically symmetric modulus 

M^y^)=v + (T)e- 2x ^ l2 e^\ (5) 

where the real function \ determines the phase of the bubble. We use this configuration 
to represent all subcritical bubbles. Subcriticality implies that the subcritical bubble 
diameter / <C R c , the critical bubble radius. Supposing that the pre-existing broken 
phase domain (bubble) is large enough compared to the subcritical bubble its wall can 
be approximated to be planar. Thus the background configuration is given by 

<S> bg (z)=v + (T)e- 2e ^ z2 / 12 , (6) 

where 9 is the step function. The formation rate of subcritical bubbles in the presence 
of $5 5 at the distance u is determined by the conditional probability and expressed as 

rV[x] ~T 4 e" /3(S3[ * 6 » + *» I " S3[ * 6 » I) , (7) 



where the superscript u denotes that the center of the bubble has been moved to the 
point x = (0, 0, u). 



Now, the normalized probability density of phase configuration x is formally given 

by 

P[x]=^- 1 T'y, (8) 

where 

and Af = JVxT' v . The evaluation of the functional integral is, of course, a hopeless 
task. In what follows we shall take some representative functions x to study statistical 
averages of some parameters using T' v as a probability density. 

To be specific, we write x = $x' where 5 is a constant angle representing the 
size of the fluctuations and x' is one °f the functions Xi( z ) — 1/2(1 + tanh(z — it)), 
X2(z) = 9(z)9(u — z)z/u + 9(z — u) or Xz( z ) = 1- Note that xi an d Xi have asymptotic 
values and 1 corresponding to the limits z — > — oo and z — > oo, respectively. The 
function X3 can be thought to be the limiting case in such a class of functions. Thus 
we can consider T' v as the probability density of the random variables u and S with 
the normalization factor 

A/ 7 = y dud8T'y[8x'} (10) 

Note, that the integrand is an even function of 5. The parameter 5 can have any values 
in the cases of xi an d Xi but by periodicity it is limited between — n and n in the case 
of X3- 

The next step is to determine the maximal distance u max . In general u is just a 
free parameter, but e.g. in the case of induced nucleation the new bubbles can not be 
formed arbitrary distant from the wall. One has to require that the new bubbles join 
to the pre-existing one forming one connected domain. Let us for the moment assume 
that the all bubbles, both the pre-existing one and the newly formed ones, have same 
phase; an assumption which certainly gives a maximal distance. Suppose also that 
the pre-existing bubble has been formed layer by layer from subcritical bubbles. By 
packing subcritical bubbles so that the minimum in the middle of the bubbles the field 
value is larger that t>_, which is required that the field at that point would roll down 
to the broken minimum, not towards unbroken one, one finds maximal possible u 
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corresponding to the tightest (lattice) packing of spheres in 3 -dimensions, face cen- 
tered cubic lattice However, at T c , where v + /v_ = 2 simple cubic lattice is equally 



good leading to same u max = 1.36/ c but for v+/v- > 2 it results a smaller one. The 
mean phase fluctuation size 5 can now be given by 



Umax 

Jd5 J du5 2 P(u,5), (12) 



and similarly for u. 

We are mostly interested in weak phase transitions where the actual transition 
temperature Tf is close to the critical temperature T c . Thus we, as an approximation, 
perform the analysis at T c except in some formulas where the tiny difference between 
T c and Tf is crucial. At small supercooling limit Tf is determined by 

a 2 9/2 ttA 3/2 

s ^ T = ^W^W^ im (13) 

leading to A = 1 — 0.0442a 1//2 /A 3//4 . Small supercooling limit is valid if 1 — A <C 1, i.e. 
a/A 3 / 2 <C 500. The minimum at Tf is given by 

a 9/2 

- e = V(v+(Tf),Tf) ~ -0.00218^^. (14) 
At the critical temperature the non-zero minimum of the potential reads 

v + (T c ) = ~T C (15) 



and the mass is given by 



m{T c f = l^-Tl (16) 



The correlation length in symmetric phase is 

lc = (17) 

and, although all sizes of subcritical bubbles exists in the symmetric phase, we use a 
representative subcritical bubble P, |10| with I = l c . 

The effect of approximation Tf ~ T c and use of I — l c is that the quantity (3 S3 
depends on the potential parameters only through the combination 7 = a/A 3 / 2 . The 
numerical evaluation of the action shows, provided that 7 is not too close to zero i.e. 
the phase transition is not too weak, that the main contribution to the averages 5 2 
and u comes from the region where S is close to zero and u is close to u max in the cases 
of functions Xi an d X2- But in the case of xs the averages receive also a significant 
contribution from the region where 5 is close to 7r and u is close to zero, as shown in 
Fig. 1. This is to be interpreted that in the case of the constant phase there are in 
average two kinds of subcritical bubbles: those having almost opposite phase as the 
pre-existing bubble and lying near the wall and those having almost the same phase 
as the pre-existing bubble and lying as far as possible from the wall. For xi an d X2 
the averages are plotted in Fig. 2. 

By comparing the average formation rate T v of subcritical bubbles near the domain 
wall to the rate T s of a single bubble not influenced by the wall given by |10[ T s oc e~ 2 ' 067 



we find that for values 7 > .5 the rate is significantly larger near the domain wall, Fig. 
3. Thus the domain wall is shielded by the subcritical bubbles nearby it. The shielding 
is even stronger for larger values of 7. This would enhance the induced nucleation 
by making subcritical bubble nucleation easier. However, the relative phases of the 
bubbles do not correlate any more and therefore we have to estimate the influence of 
the phase differences. 

The influence of phase deviations of the size 5 to the bubble dynamics can be most 
easily analyzed within the context of thin wall approximation valid in weak phase 
transition case. Thin wall approximation requires that the maximum of the potential 
Vmax 3> e, which implies a/ A 3 / 2 <C 10 in accordance with small supercooling limit. 
The critical bubble radius R is obtained by extremizing the bounce action 



B = --ne eff R 3 + AnaR 3 , (18) 
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where the surface tension in thin wall approximation is given by 



o = (19) 



Usually e e ff is simply given by Eq. fll4j) but now it is affected by phase fluctuations 
in the broken phase, so that 

e eff = e - Ae, (20) 

where 

^R 3 Ae=( J d 3 x\V<S> R (x)\ 2 ). (21) 

xeB(Q,R) 

is now a random configuration having modulus v + and a varying phase between 
values —5 and S. 

If we suppose that the typical size of the spatial extension of the phase fluctu- 
ations is u, as can be done in the cases of Xi an d X2, we are able to calculate 
Ae. When the maximum values of phase fluctuations are situated at the points 
x = (n x u,n y u,n z u), n x , n y , n z e Z the calculation can be reduced to a unit lat- 
tice where n x , n y , n z = 0, 1. Taking = e tSx ^ the expectation value Ae can be 
now expressed in terms of the unit lattice: 

lrfAe = ^ 2 (|) 3 ( / rf 3 x(Vx) 2 ). (22) 



xG[0, u] 
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At each corner we associate equally distributed independent random variables fn x ,n y ,n z 
with values ±1 and correlations (fnnn) = ■ Moreover, we choose the phase 



configuration in the unit bravais lattice to be the simplest possible polynomial, i.e. 
the Lagrangian interpolation polynomial 



n x ,n y ,n z =0,l k x ,k y ,k z =0,1 " " 

kx 7^ Tlx 
k y 7= n y 
fc 2 / n 2 



An easy calculation yields now 



Ae = ^±P. (23) 
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Because from Eq. ( [TED one obtains the critical bubble radius 

i? c = — (24) 

€ eff 

the critical bubble radius tends to infinity if e e ff — > 0. Hence a general requirement 
for existence of the critical bubble can be stated to be e e ff > 0, i.e. 

*< 0.129- ( w ) . (25, 

Thus Eq. (p5]) expresses an overall upper bound for mean phase fluctuations inside 
a critical bubble during the cosmological electroweak phase transition in terms of 
phenomenological parameters and fluctuation correlation length. Comparing Eq. (p5|) 
and the Fig. 2 one can conclude that it is not possible for any size of critical bubble 
to survive from the effects of phase fluctuations. One should have a mechanism which 
would damp more efficiently the variations of 5. 

In the case of X3 the energy change has to be calculated differently. Because the 
small, opposite phase 6 = it bubbles can be approximated to be isolated in the 5 = 
phase, the energy change is approximated by Se = , where E2 is the energy of 
single wrong-phase bubble and V\ and t>2 are the volumes of 5 = and 5 = tc bubbles, 
correspondingly. A calculation yields an lower bound for 5e > 0.05957. Requiring 
e e fj > we obtain aX 1 ^ 2 > 0.55 or, because 7 10, A 3> 0.23 which can hardly be 
true. Hence also %3 case confirms the conclusion that the phase deviations are too 
large for induced nucleation to be possible. 

In the present paper we have analyzed the subcritical bubble nucleation near a 
domain wall. The wall itself induces nucleation nearby it so that the subcritical bubble 
formation rate is larger close the wall than far of it. It also appears that deviations of 
the phase of the field can be remarkably large and prohibits the system from induced 
nucleation. Note, that if the size of typical subcritical bubble is smaller that the 
correlation length as suggested lately |L2] , the bound (|25| ) would be even stricter. Also 



those subcritical bubbles having very small amplitude compared to v + may have effect 



to the large amplitude configurations. In realistic theories, like the Standard model 
where the symmetry group is larger, the actual numbers for 5 and u would surely alter 
but, because they have more degrees of freedom, the qualitative results are likely to 
hold. 
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Figure captions 



Figure 1. The fraction of 500 randomly generated points in the 8-u plane which pass 
the limit @S 3 < 1 with 7 = 1. Open circles correspond to the case of xi an d filled 
circles to the case of X3- 

Figure 2. The averages S[xi] (solid line), 5[x2\ (dashed line), u[xi] (dotted line) and 
u[x2\ (dash-dotted line) as a function of 7. For 7 = the average 8 diverges as a 
function of the upper limit of the (^-integration. The upper limit used in the plot is 
8 m ax — 10 2 . Studying the dependence of the upper limit of the integration numerically 
reveals that it is wery weak for 7 > 0.1, implying that only those values of 5 with 7 
larger than that limit are reliable. Note that 7 = implies that the scalar potential 
has only the symmetric phase minimum |$| =0. 

Figure 3. The average formation rate T' v of sub critical bubbles near the domain wall 
(solid line) and the rate T s = exp(— 2.O67) of a single bubble not influenced by the 
wall (dashed line). 
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